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LAMINATION RESIDUAL STRESSES IN MULTILAYERED FIBER COMPOSITES 


by Christos C. Chamis 
Lewis Research Center 

SUMMARY 

Residual stresses arising from the lamination fabrication process of composites 
are investigated using linear laminate theory. The resulting equation for predicting 
residual stresses is presented. The variables that influence the residual stresses are 
identified. These variables are constituent materials and their thermal and elastic 
properties, average fiber diameter and number of fibers per end, fiber and void vol- 
ume ratios; ply stacking sequence and orientation, and processing temperature. The 
effects of these variables on the residual stresses of several composite systems are 
examined. The fiber composites considered are boron/aluminum, boron/epoxy, 

S- glass/epoxy, Thornel-50 graphite/epoxy, and Modmor-I graphite/polyimide. 

It was found that the residual transverse tensile and in- plane shear stresses can 
reach magnitudes comparable to the corresponding ply strengths and can cause transply 
cracking. The ply residual stresses are very sensitive to constituent material proper- 
ties, ply stacking sequence and orientation, fiber content, and processing temperature. 
They are relatively insensitive to the void content. The residual stresses are relatively 
small in composites where the fiber content or the processing temperature vary only 
moderately from ply to ply. Ply residual stresses increase with increasing matrix 
modulus. They increase linearly with increasing matrix thermal coefficient of expan- 
sion. It is possible to orient plies in simulated components so that some plies have 
tensile longitudinal residual stress. Residual stresses can cause interply delamination 
due to the relative rotation of adjacent plies. 

Residual stress transply cracking, in general, can be prevented by increasing the 
ply transverse tensile and shear strengths, decreasing the processing temperature, or 
decreasing the difference between the fiber and matrix thermal coefficients of expansion. 
In specific designs transply cracking can be prevented by suitable ply orientation, in- 
crease in fiber content, use of transitional plies, and reduction of the number of off- 
axis plies. 

Experimental work to fully assess the residual stress effect on the structural re- 
sponse of the composite is suggested. 



INTRODUCTION 


Multilayered fiber composites are readily adaptable to simultaneous component and 
material designs. This approach allows the designer to take full advantage of the high 
stiffness and strength-to-density ratios of fiber composites. The result is an efficient 
use of fiber composite material in structural components. The composite is usually 
designed by specifying the fiber direction in each ply, the number of plies, and the fiber 
volume ratio necessary to meet the given loading requirements. 

However, having different fiber directions in different plies can cause transply 
cracks such as have been observed in graphite -fiber composites. These cracks are 
caused by the large thermal residual stresses present in the composite. The residual 
stresses result from the fabrication process used to make fiber composites. The re- 
sidual stresses arising from the fabrication process are examined quantitatively herein 
using linear laminate theory. 

Multilayered fiber composites consist of several unidirectional layers (plies) 
oriented at specified angles relative to some coordinate axes system. The geometry of 
such a composite is depicted in figure 1. The plies consist of parallel fibers embedded 
in a matrix. The geometry of a ply is shown in figure 2. Glass or graphite fiber com- 
posites will have many filaments through the ply thickness. 

In general, composites achieve their structural integrity as follows: Nonmetallic 
matrix composites are cured at elevated temperature under pressure; metallic matrix 
composites are either hot pressed or plasma sprayed. These methods induce residual 
stresses in the composite when it is at a temperature different from its cure or proc- 
essing temperature. The difference between the cure or processing and the use tem- 
perature will be referred to here as temperature difference. The residual stresses can 
be either microresidual or macroresidual stresses. 

The microresidual stresses are influenced by the local geometry of the constituents 
within the ply. They arise from the difference in the thermal coefficients of expansion 
of the constituents and the temperature difference. The matrix phase change shrinkage 
for resins also contributes to the microresidual stress. Adams, Doner, and Thomas 
(ref. 1) and Haener, Ashbaugh, Chia, and Feng (ref. 2), for example, treat this type of 
residual stress analytically. Marloff and Daniel (ref 3) and Koufopoulos and Theocaris 
(ref. 4) examined microresidual stress experimentally. The ply shown in figure 2 will 
have microresidual stresses. 

Macroresidual stresses (ref. 5) are integrated averages through the ply thickness. 
These stresses have an effect on the ply similar to an equivalent thermal or mechanical 
load applied to the composite. They arise from the difference between the ply longitu- 
dinal and transverse thermal coefficients of expansion. Macroresidual stresses are 
present in those composites that have plies oriented at different angles and whose use 
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temperature differs from the cure or processing temperature. These stresses will be 
present in the composite depicted in figure 1. 

The significance of the lamination residual stress was realized when transply 
cracks were observed in 0°/90° graphite composites. Doner and Novak (ref. 5) pre- 
sented photomicrographs showing transply cracks. DeCrescente and Novak (ref. 6) and 
Winters (ref. 7) have observed similar cracks in other graphite/epoxy composites. 

Some analytical treatment of the lamination residual stress has also been reported. 
Doner and Novak (ref. 5) considered the residual stresses in 0°/90° composites. 
Edighoffer, Ravenhall, and Juneau (ref. 8) examined the residual stresses in a boron/ 
epoxy composite with laminations ± 9 / 90 where the ± 6 are considered as an equivalent 
0° ply. This author (ref. 9) used linear laminate theory to predict residual stresses in 
graphite composites without restriction to ply orientation. The lamination theory used 
in reference 9 was developed in reference 10. It was implemented in references 11 
to 13 and in a computer code (ref. 14). 

The computer code of reference 14 provides a method for micromechanics, lami- 
nation, and stress analyses of multilayered fiber composites. The inputs to this code 
are constituent material properties, fiber and void contents, correlation coefficients, 
composite geometry, temperature difference (use temperature minus cure tempera- 
ture), and composite stress resultants or displacements. 

The computer code of reference 14 permits the variation of several basic param- 
eters in fiber composites. It is used herein to investigate the residual stresses in lami- 
nated fiber composites. The following systems are examined: boron/aluminum, boron/ 
epoxy, S-glass/epoxy, Thornel-50 graphite/epoxy, and Modmor-I graphite/polyimide. 
The significant parameters in these systems are ply orientation, fiber and void contents, 
matrix modulus of elasticity, matrix and thermal coefficients of expansion, and the tem- 
perature difference. To the author T s knowledge, the effects of some of these param- 
eters on the residual stresses are investigated for the first time. Also investigated for 
the first time are (1) possible ply de lamination due to residual stress, and (2) residual 
stresses in composites which are unsymmetric with respect to bending. 

The ply orientations examined are shown in table I together with the notation adopted 
herein for convenience. Theoretical thermal and elastic properties for the composite 
systems considered are listed in table II. These properties are for a fiber volume ratio 
of 0.55. The corresponding theoretical strength properties are given in table III. The 
results in tables n and III are those predicted by the computer code (ref. 14). 

The significant parameters mentioned previously are varied in the computations. 

The results are plotted to illustrate the dependence of the ply residual stresses on these 
parameters as is described subsequently. The plots are grouped by composite system. 
The composite systems boron/aluminum, boron/epoxy, S-glass/epoxy, and Thornel-50 
graphite/epoxy are used in some basic ply orientations. Results for various ply orien- 
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tations in composites unsymmetric with respect to bending are also presented for these 
systems. The composite system Modmor-I graphite/polyimide is used in simulated 
structural components. Some results are presented for this system to demonstrate 
possible ply de lamination due to residual stresses. 

The results presented herein are for three temperature differences: 900° F (773 K) 
for boron/aluminum, 300° F (440 K) for the boron/epoxy, S- glass/epoxy, and 
Thornel-50 graphite/epoxy systems, and 600° F (606 K) for the Modmor-I graphite/ 
polyimide system. Residual stresses for any other temperature differences are ob- 
tained from the plots by direct proportion. This is possible because a linear theory is 
used and the residual stress is linearly dependent on the temperature difference. 

Several recommendations are made for minimizing or eliminating the residual 
stress. These recommendations are deduced from the results presented and the equa- 
tions used to predict the residual stresses. These equations are presented in their 
final form. The derivations are given in reference 10. 


SYMBOLS 

A array of composite extensional stiffnesses referred to composite structural 

La 

axes 

C cx array of composite coupling stiffnesses referred to composite structural axes 
D cx array of composite bending stiffnesses referred to composite structural axes 

in 

array of strain-stress relations (elastic constants of the i 1 ply) 

E nl ply longitudinal modulus 

E ^2 P^y transverse modulus 

G^i 2 ply (in-plane) shear modulus 

kj,k y fiber and void volume ratios, respectively 

^cATx vec t° r of unbalanced thermal moments referred to composite structural axes 

N cATx vec t° r of unbalanced thermal forces referred to composite structural axes 
array of transformation coefficients for the i Ln ply 
S^ht ply longitudinal tensile strength 

S^nc ply longitudinal compressive strength 

£>Z22T ply transverse tensile strength 

£>Z22C ply transverse compressive strength 
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S Z12S 

AT Zi 

x,y,z 

z Zi 

1,2,3 

a Zi 

“Zll 

a Z22 

%2T 

e cox 


e Zi 

e mpT 

6 


K cx 
v l 12 


CT Zi 

^M22 
{ } 

[ ] 


ply in-plane shear strength 
temperature difference for the i ply 
structural axes 

XT* 

distance from reference plane to centroid of the i ply 
material axes 

XT. 

vector of thermal coefficients of expansion of the i u ply referred to ply's 
material axes 

ply longitudinal thermal coefficient of expansion 
ply transverse thermal coefficient of expansion 
correlation coefficient for ply transverse tensile strength 
void strain magnification coefficient 

vector of composite strains referred to composite structural axes at the ref- 
erence plane 

XT* 

vector of strains for the i Ln ply referred to ply's material axes 
matrix allowable strain in tension 

ply angle measured from the composite structural axes to the ply material 
axis 1 

composite local curvatures referred to composite structural axes 

ply Poisson's ratio in the two direction when the load is applied in the one 
direction 

vector of ply stresses referred to ply's material axes 
strain magnification factor for ply transverse strength 
vector or column matrix 
square array 


GOVERNING EQUATIONS 

The laminate theory equation for predicting the ply residual stresses is given by 

{cr a } = [E^] 1 <[R a ] {e cox } " z Zi[ R Zi^ K cx* -AT Zi^ Q! Zi^ (1) 

The reference plane strains {e } and the curvature changes { k } for a free com- 

v UA LA 
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posite (no external loads and no boundary constraints) are computed from 


pw] 

f A cx][ C cx] 

-1 

KaTxO 

1 r 


< 

! 

j 

KcU 

[ C cx][ D cxi 


KatxU 


( 2 ) 


Equations (1) and (2) show that the ply residual stress is a function of the following 
factors: (1) the composite structural stiffness, (2) the ply spatial location in the com- 
posite, (3) the ply stress-strain relations, (4) the ply thermal coefficients of expansion, 
and (5) the temperature difference between the ply and the reference value. This differ- 
ence equals the ply temperature minus the cure or processing temperature in computing 
residual stresses. The factors (1), (3), and (4) are related to the basic variables (con- 
stituent materials and thermoelastic properties, fiber and void volume ratios, number 
of plies, and ply orientation). These relations are given in either reference 10 or 14. 


RESIDUAL STRESSES IN BORON/ALUMINUM COMPOSITES 

The temperature difference selected for boron/aluminum composites is 900° F 
(773 K). The room temperature properties are used in the linear analysis. The appli- 
cability of linear analysis and room-temperature properties might be questioned for this 
temperature difference, particularly for the aluminum matrix. However, there are four 
reasons that justify this approach. First, the residual stresses are computed at room 
temperature. Second, the aluminum matrix modulus decreases with increasing temper- 
ature, but its thermal coefficient of expansion increases. These two effects tend to be 
compensatory. The third reason is that the results from a linear analysis can be scaled 
down proportionally to any desired temperature difference. The fourth reason is that 
composites that are designed to carry the residual stress at room temperature caused 
by a 900° F (773 K) temperature difference will be conservatively designed. 

Ply residual transverse and shear stresses (as defined in fig. 2) are plotted as 
functions of ply orientation angle and fiber volume ratio in figure 3 . The composite 
geometry and the temperature difference is noted in the captions of these figures. 

Ply residual transverse stresses as a function of ply angle are plotted in figure 4 
for the general composite geometry 8[0, 2(±0), 2(90), 2(+0), 0]. The temperature dif- 
ference and the fiber volume ratio (FVR) are 900° F (773 K) and 0. 55, respectively. 

The results in this figure for 6 - 0° or 90° illustrate the dependence of the ply residual 
stress on the number of plies in each orientation. 

Ply residual transverse stresses for composites of the following three geometries 
8[4(+0), 4(-0)], 8(±0), 8[4(±0), 4 (t0)] and with an FVR of 0. 55 are plotted against 0 in 
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figure 5. The ply considered in these composites is the first ply in the sequence. The 
composite geometry 8[4(+0), 4(-0)] represents a composite that is nonsymmetric with 
respect to bending. The composite geometry 8(±0) represents a composite that is non- 
symmetric with respect to bending with interspersed plies. The composite geometry 
8[4(±0), 4(+0)] represents a composite that is symmetric with respect to bending with 
interspersed plies. 

The results in figures 3 to 5 show that the ply residual transverse and shear 
stresses increase with ply angle but decrease with fiber volume ratio. The smaller the 
number of plies in a given orientation, relative to other orientations in the same com- 
posite, the greater the residual stresses in these plies. The residual stresses in the 
outer plies of nonsymmetric noninterspersed composites are quite high. They are about 
twice as high as those in either nonsymmetric interspersed or symmetric interspersed 
plies. Introduction of ±45° plies for increased torsional stiffness produces high residual 
stresses. The residual stresses in quasi (pseudo)-isotropic symmetric composites (such 
as 8[0, 2(±45), 2(90), 2 (+45), 0], fig. 6) are equal in all the plies. 

The transverse ply tensile strength for a boron/aluminum unidirectional composite 
(ply) at an FVR of 0. 55 is 14. 1 ksi (97. 1 MN/m^) (see table III). This value is exceeded 
in many of the composites shown in figures 3 to 5. 

To what extent the residual stresses can be annealed, if any, and how much these 
stresses affect the in-situ ply strength is not known. Also, the effect of creep in re- 
lieving residual stresses is not known at this time. 


RESIDUAL STRESSES IN BORON, S- GLASS, AND 
THORNEL-50 GRAPHITE /EPOXY COMPOSITES 

The residual stresses in these composites are computed from a 300° F (440 K) tem- 
perature difference, room temperature constituent material properties, and the linear 
theory available in reference 14. This theory is briefly described in the INTRODUC- 
TION. Five reasons justify this approach: (1) The residual stresses are computed at 
room temperature. (2) Polymerization shrinkage (matrix phase change), though not 
accounted for directly in the theory, is compensated for by the higher temperature. 

(3) The epoxy modulus decreases with increasing temperature, but the epoxy thermal 
coefficient of expansion increases. (4) The residual stress results presented, herein, 
can be scaled down proportionally to any desired temperature difference. (5) Com- 
posites designed to withstand the residual stress at room temperature predicted by this 
approach will be conservatively designed. 

Results for ply residual transverse stresses are presented. These stresses are 
primarily responsible for the transply cracking. 
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Boron/Epoxy Composites 


The residual transverse stresses in the ±0 plies of the boron/epoxy composites 
8[2(0), 2 (±0), 2(+0), 2(0)] are plotted against 0 and three fiber volume ratios in fig- 
ure 6(a). 

Ply residual transverse stresses for the 0°, ±0°, and 90° plies of the composites 
8[0, 2(±0), 2(90), 2 (+&), 0 ] are plotted against 0 in figure 60b). The FVR is 0. 55 for these 
composites. The results in this figure at 0 = 0° and 90° illustrate the effect on re- 
sidual stress of the numbers of plies in each orientation. 

The residual transverse stress in the outer ply of the composites 8[4(+0), 4(-0)], 
8(±0), and 8[4(±0), 4(+0)] is plotted against 0 in figure 6(c). The FVR is 0. 55 for these 
composites. The composites 8[4(+0), 4(- 0)] are nonsymmetric with respect to bending 
and have a noninterspersed ply stacking sequence. The composites 8{±0) are nonsym- 
metric with respect to bending but do have an interspersed ply stacking sequence. The 
composites 8[4(±0), 4(+0)] are symmetric and have interspersed ply sequence. 

The transverse strength for the unidirectional composite (ply) at an FVR of 0. 55 is 
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8. 6 ksi (59. 3 MN/m ) (see table III). This value is exceeded only by the residual 
stresses in the outer plies of the composite 8[4(+0), 4(-0)] (fig. 6(c)). However, the 
margin of safety is about 0. 3 in many of the other composites. 

Several interesting points are illustrated by the results shown in figure 6: (1) The 
ply residual transverse stress in boron/epoxy composites is relatively high when com- 
pared with the ply transverse strength for some composite geometries. (2) The ply re- 
sidual stress is quite sensitive to the ply angle, but not as sensitive to FVR. (3) The 
ply transverse residual stresses in the composites 8[0, 2(±0), 2(90), 2(+0),O] are practi- 
cally invariant with 0. (4) Nonsymmetric noninterspersed- ply composites experience 

high ply residual stresses when compared with other composite geometries. (5) Non- 
symmetric interspersed- ply and symmetric interspersed- ply composites experience ply 
residual stresses of about the same magnitude. 


S-Glass/Epoxy Composites 

The residual transverse stress in the ±0 ply of the S-glass/epoxy composites 
8[2(0), 2(±0), 2(+0), 2(0)] is plotted against 0 and three FVR in figure 7(a). 

Ply residual transverse stresses for the 0°, ±0°, and 90° plies of the composites 
8[0, 2(±0), 2(90), 2(+0), 0] are plotted against 0 in figure 7 (b). The FVR is 0. 55 for 
these composites. The results in this figure for 0-0° and 90° illustrate the effect on 
residual stress of the number of plies in each orientation. 

The residual transverse stresses in the outer ply of the composites 8[4(+0), 4(-0)], 
8(±0), and 8[4(±0), 4 (+0)] are plotted against 0 in figure 7(c). The FVR is 0. 55 for 
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these composites. The composites 8[4(+0), 4(-0)] are nonsymmetrie noninterspersed- 
ply composites. The composites 8(±0) are nonsymmetrie interspersed-ply composites. 
The composites 8[4(±0), 4(+0)] are symmetric interspersed-ply composites. 

The transverse strength for the S-glass/epoxy unidirectional composite (ply) at an 
FVR of 0. 55 is 4. 4 ksi (30. 4 MN/m^) (see table III). This value is very close to the ply 
residual stress of several composites. It is exceeded in the outer plies of some com- 
posites 8[4(+0), 4(-0)]. This implies that S-glass/epoxy composites of the geometries 
and FVR examined in figure 7 could experience transply cracks . This will not be the 
case if the transverse strength of the plies is much greater than 4. 4 ksi (see, e. g. , 
table I of ref. 5). 

Several interesting points are illustrated by the results in figure 7. (1) The calcu- 

lated ply residual transverse stress could exceed the ply transverse strength. The ex- 
tent that the ply transverse strength is exceeded by the calculated ply residual stress 
depends on composite geometry. (2) The ply residual stress is quite sensitive to ply 
orientation angle in certain composite geometries. (3) The ply residual stress is rela- 
tively insensitive to fiber volume ratio in composites with less than 0. 55 FVR. (4) The 
ply residual stress is sensitive to the number of plies in each orientation. (5) Nonsym- 
metric noninterspersed-ply composites experience comparatively large ply residual 
stresses. (6) Nonsymmetrie interspersed-ply and symmetric interspersed-ply com- 
posites experience ply residual stresses of about the same magnitude. 


Thornel-50 Graphite/Epoxy Composites 

The residual transverse stress in the ±0 ply of the Thornel-50 graphite/epoxy 
composites 8[2(0), 2(±0), 2(+0), 2(0)] is plotted against 0 and three values of FVR in fig- 
ure 8(a). 

Ply residual transverse stresses for the 0°, ±0°, and 90°-plies of the composites 
8[0, 2(±0), 2(90), 2(+0), 0] are plotted against 0 in figure 8(b). The FVR is 0. 55 for these 
composites. The results in this figure for 0=0° and 90° illustrate the effect on re- 
sidual stress of the number of plies in each orientation. 

The residual transverse stress in the outer ply of the composites 8[4(+0), 4[-0)], 
8(±0), and 8[4(±0), 4(+0)] is plotted against 0 in figure 8(c). The FVR is 0.55 for these 
composites. The composites 8[4(+0), 4(-0)] are nonsymmetrie noninterspersed-ply 
composites. The composites 8(±0) are nonsymmetrie interspersed-ply composites. 

The composites 8[4(±0), 4(+0)] are symmetric interspersed-ply composites. 

The transverse strength for the Thornel-50 graphite /epoxy unidirectional composite 
(ply) at an FVR of 0. 55 is 3. 7 ksi (25. 5 MN/m^) (see table III). This value is exceeded, 
or nearly so, by several composites. This implies that plies in Thornel-50 graphite/ 
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epoxy composites of the geometries and values of FVR examined in figures 8(a) to 8(c) 
will experience transply cracks. This has been observed experimentally (refs. 5 to 7). 
This will not be the case if the transverse strength of the plies is about 8. 0 ksi 
(55. 2 MN/m 2 ). 

Several interesting points are illustrated by the results in figure 8. (1) The calcu- 

lated ply residual transverse stress exceeds the ply transverse strength. The extent 
that the ply transverse strength is exceeded by the calculated ply residual stress de- 
pends on composite geometry. (2) The ply residual stress is quite sensitive to ply angle 
in certain composite geometries. (3) The ply residual stress is sensitive to FVR. 

(4) The ply residual stress is insensitive to the number of plies in each orientation. It 
is relatively insensitive to 6 in the composites 8[0, 2(±0), 2(90), 2(+0), 0]. (5) Nonsym- 

metric noninterspersed-ply composites experience ply residual stress that is about 
30 percent higher than other composite geometries. (6) Nonsymmetric interspersed- ply 
and symmetric interspersed-ply composites experience ply residual stresses of about 
the same magnitude. 


SIMULATED COMPOSITE COMPONENTS 

The residual stresses in the plies of two simulated composite components were 
studied. The composite geometries are 24[4(±30), 2(±15), 12(0), 2(+15), 4(+30)] and 
24[4(±45), 2(±22. 5), 12(0), 2(+22. 5), 4(+45)]. The composite system is Modmor-I graphite/ 
polyimide, and it is used for elevated temperature application. The temperature differ- 
ence for this system is about 600° F (606 K). 

The dependence of the ply residual stresses on several factors is examined. These 
factors are (1) fiber volume ratio (FVR), (2) void content, (3) variation of the FVR from 
ply to ply within the same composite, (4) different groups of plies with different FVR, 

(5) variations in the matrix thermal coefficient of expansion, and (6) variations in the 
matrix modulus. Also, the effect of the residual stresses on the relative rotation of 
adjacent plies is investigated. 

The computations are based on a 600° F (606 K) temperature difference and room- 
temperature constituent material properties. The justification for this approach is the 
same as that given for the other fiber/epoxy composite systems. (See section Residual 
Stresses in Boron, S-Glass, and Thornel-50/Epoxy Composites. ) 

The ply residual transverse stress is plotted against FVR in figure 9(a) for both 
components. Results for the ply residual longitudinal stress are shown in figure 9(b). 

The residual longitudinal stress in the 0° plies of both components is tensile. Results 
for ply residual shear stress are plotted against FVR in figure 9(c). The results in 
these figures show that the ply residual stresses are sensitive to both composite geom- 
etry and FVR. 
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The void content effect on the ply transverse residual stress is illustrated in fig- 
ure 10. This effect is negligible for all practical purposes. However, it is important 
to keep in mind that the ply transverse strength decreases with increasing void content 
(ref. 12). The decrease in ply transverse strength could produce transply cracks at 
smaller temperature differences. Proper consideration of the void effects is necessary 
to avoid transply cracking. (See eq. (3) in this section. ) 

The effects of variable ply FVR on the ply transverse residual stresses are shown 
in figure 11 for one component only. These effects are negligible as can be seen in the 
figure. Having different groups of plies at different FVR has a negligible effect on the 
ply residual transverse stress. This is illustrated in figure 12. Negligible effects were 
found also when the ply temperature was varied similar to the FVR. These results are 
not presented here. 

The dependence of the ply residual transverse stress on the matrix thermal coeffi- 
cient of expansion is illustrated in figure 13. Two important points are observed here: 

(1) The ply residual stress varies linearly with matrix thermal coefficient of expansion 
and (2) the residual stresses can be reduced by decreasing the magnitude of the thermal 
coefficient of the matrix. The effects of varying the fiber thermal coefficients of expan- 
sion will produce an effect opposite to those shown in figure 13 for the thermal coeffi- 
cient of the matrix. This implies that increasing the fiber thermal coefficient of expan- 
sion is another effective way to reduce the residual stress. 

The dependence of the ply residual transverse stress on the matrix modulus is illus- 
trated in figure 14. The results in this figure show that the residual stress increases 

nonlinear ly with increasing matrix modulus. This increase is quite rapid at modulus 

o 

values of 100 to 700 ksi (690 MN/m ) and not as rapid at modulus values greater than 
700 ksi (4. 8 MN/m 2 ). 

The results in figure 14 should be examined jointly with the ply transverse strength 
before any meaningful conclusions are made. The ply transverse strength is given by 
the expression (ref. 12). 


S i22T = %2T E m (3) 

PyV ju22 

where S^2 t is the transverse tensile strength, @22T * s the corre l a tion factor, 
e mpT is the a ll° wa kle matrix tensile strain to denote failure, 0 V is the void strain 
magnification factor, <P ^22 is the s t raan magnification factor, and E ^2 the ply 
transverse modulus. The variables e m p>p> <P^ 22 ’ anc * ^“122 are not independent of the 
matrix modulus. For example, <P ^22 decreases but E ^2 increases with increasing 
matrix modulus (ref 10), and e m p T increases with increasing matrix modulus in gen- 
eral. 
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The following point can be made in view of the preceding discussion. The ply trans- 
verse residual stresses increase with increasing matrix modulus. The ply transverse 
strength also increases with increasing matrix modulus. Transply cracking will even- 
tually occur if the ply residual transverse stress increases at a higher rate than the ply 
transverse strength. Limited results showed that the ply transverse strength increases 
at a higher rate than the residual stress for some graphite composites during this inves- 
tigation. 

The effect of the residual stress on the adjacent ply relative rotation is illustrated 
in figure 15. Relative rotation is plotted against fiber volume ratio for selected adjacent 
plies of the two components. The estimated allowable relative rotation is also plotted in 
the figure, and the prediction of it is discussed in references 10 and 11. The results in 
figure 15 show that residual stresses could cause interply de laminations. The de lami- 
nation depends on the constituent materials and composite geometry, the ply angles of 
the two adjacent plies, and the fiber and void volume ratios. It is interesting to note 
that the component with the ±30° and ±15° plies has a higher probability of interply de- 
lamination than the component with the ±45° and ±22. 5° plies. The component with the 
±30° and ±15° plies is superior in all other respects with respect to residual stress. To 
the author's knowledge, the significance of the results in figure 25 has not been recog- 
nized by other investigators to date. 

The results in figures 16 to 25 show that composites designed by neglecting the re- 
sidual stresses will be unconservative in some loading conditions. 

The effects of introducing transitional plies (intermediate plies for less-abrupt 
orientation transition between groups of plies) or decreasing the number of off-axis plies 
on the residual stress are illustrated in figure 16. The residual transverse stress in 
the ±45° plies is plotted against FVR. The composite with the transitional plies is de- 
noted by the long-dashed curve. The composite in which the transitional plies have been 
incorporated with the face (surface) plies is denoted by the solid curve. The composite 
in which the transitional plies have been incorporated in the core plies is identified by 
the short- dashed curve. 

The results in figure 16 show that both reduction of the number of off-axis plies and 
introduction of transitional plies reduce the residual stress by small amounts. The 
number of off-axis plies is controlled by strength and stiffness requirements and changes 
might not be permitted in specific designs. 


WAYS TO MINIMIZE OR ELIMINATE LAMINATION RESIDUAL STRESSES 

Several alternatives are possible for minimizing or eliminating the lamination re- 
sidual stresses. These alternatives may be grouped into two categories. The first 
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category is general, and the second is for specific designs. They are summarized 
briefly. 

In general, residual stresses at room temperature can be reduced or eliminated by 

(1) Use of lower cure or processing temperature matrices 

(2) Decreasing the difference between the fiber and the matrix thermal coefficients 

of expansion 

(3) Increasing the allowable elongation of the matrix at failure with a simultaneous 

increase in the matrix modulus. 

In specific designs, the residual stress can be reduced by 

(1) Proper selection of composite system 

(2) Suitable ply orientation and stacking sequence 

(3) Appropriate selection of the number of plies in each direction 

(4) Introduction of transitional plies 

(5) Decreasing the number of off-axis plies if stiffness and strength requirements 

permit it 

(6) Increasing the fiber volume ratio 

(7) Any possible combination of these. 


EXPERIMENTAL WORK NEEDED 

Both theory and experiment show that the ply residual transverse stresses due to 
lamination are of magnitudes comparable to those of the ply transverse strengths. How- 
ever, four important questions need to be answered before the effect of residual 
stresses on the composite structural response can be realistically assessed: (1) decay 
of residual stress with time, (2) effect of transply cracks on composite stiffness, 

(3) effect of transply cracks on composite strengths, and (4) effect of transply cracks on 
composite degradation by the environment. 

These questions can be answered best by simple but strategic experimental investi- 
gation. Bending tests in fatigue and creep should provide qualitative answers to these 
questions. 


CONCLUSIONS 

The results of this theoretical investigation, which is based on linear laminate 
theory and on a semiempirical micromechanics theory, lead to the following conclu- 
sions. 

1. The ply residual transverse stress is generally tensile. The residual stress 
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depends on the composite system, ply orientation, fiber volume ratio, and processing 
conditions. The residual stress can be large relative to transverse strength for 
graphite/epoxy and boron/aluminum composites. 

2. The ply residual in-plane shear stress can be comparable to the ply's shear 
strength. 

3. The ply residual transverse or shear stress can be comparable to the corre- 
sponding ply strengths in simulated components. This can cause transply cracking. It 
can also reduce the ply capacity to carry transverse tensile or in- plane shear loads. 

4. The ply residual stresses are very sensitive to temperature difference, ply 
orientation, number of plies in each direction and to fiber volume ratio. 

5. The ply residual stresses are relatively insensitive to void content and to mod- 
erate ply- to- ply fiber volume ratio and processing temperature variations. 

6. The ply residual stresses generally increase with increasing ply orientation 
angle. 

7. The ply residual stresses increase nonlinear ly with increasing matrix modulus. 
They increase linearly with the matrix coefficients of expansion. 

8. It is possible to orient plies in simulated components so that some plies have 
tensile longitudinal residual stress. However, the residual longitudinal stress is 
usually compressive. Its magnitude is usually 15 to 20 percent of the ply's correspond- 
ing compressive strength. 

9. Transply residual stress cracking, in general, can be prevented by increasing 
the ply transverse and shear strengths, and decreasing the processing temperature. It 
can be prevented by decreasing the difference between the thermal coefficients of the 
fiber and the matrix sufficiently. 

10. In specific designs, cracking can be prevented by suitable ply orientation, var- 
iation of the number of plies in each direction, introduction of transitional plies, reduc- 
tion of the number of off-axis plies, increasing the fiber content (within practical limits), 
or by any combination of these . 

11. Neglecting the residual stresses in designing with fiber composites can lead to 
unconservative design. 

12. Experimental work is needed for a realistic assessment of residual stresses on 
composite structural response. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 5, 1970, 

129-03. 
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TABLE I. - COMPOSITE GEOMETRY AND NOTATION IDENTIFICATION 


Composite 

Ply number 

Notation 

1 

2 

3 

4 

5 

6 

7 

8 


Ply orientation angle, deg 


1 

0° 

0° 

+ 6 

-9 

-0 

+0 

0° 

0° 

8[2(0), 2(±0), 2(+0), 2(0)] 

2 

0° 

+ 9 

-9 

90° 

90° 

-0 

+ 0 

0° 

8[0, 2, (±0)2(90), 2(T0), 0] 

3 

+ 0 

-9 

+9 

-0 

-0 

+0 

-0 

+0 

8[4(±0), 4(+0)] 

4 

+ 0 

-9 

+9 

-9 

+ 0 

-0 

+ 0 

-9 

8(±0) 

5 

+ 6 

+9 

+9 

+9 

-9 

-0 

-0 

-0 

8[4(+0), 4(-0)] 
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TABLE II. - THEORETICAL UNIAXIAL PLY ELASTIC AND THERMAL PROPERTIES 


OF VARIOUS COMPOSITES AT ROOM TEMPERATURE 
(Normalized at 0. 55 fiber volume ratio and zero voids (refs. 12 and 13). ) 


Composite system 


Elastic modulus 


Poisson’s 

Thermal coefficient of expansion 


Longitudinal 

Transverse 

Shear 

ratio 

Longitudinal Transverse 



psi 

GN/m 2 

psi 

GN/m 2 

psi 

GN/m 2 


in. /in. /°F 

cm/cm/K 

in./in./°F 

cm/cm/K 

Boron/aluminum 

37.5x10® 

258. 6 

26. 9' 10 6 

185.5 

13. 2x10® 

91.0 

0.23 

4.0x10"® 

7.2X10" 6 

6.2X10' 6 

11.2x10"® 

Boron/epoxy 

33. 2x10® 

228. 9 

2. 1-10® 

14.5 

1.0x10® 

6.9 

.24 

3.0x10"® 

5.4x10"® 

14.2x10"® 

25.6x10"® 

S-glass/epoxy 

7. 0x10® 

48.3 

1.7v 10® 

11.7 

. 9x10® 

6.2 

.25 

3.7x10"® 

6.7x10"® 

14. 1x10"® 

25. 4x10"® 

Thornel-50 graphite /epoxy 

27. 8x10® 

191.7 

1. 0-- 10® 

6.9 

.7x10® 

4.8 

.24 

-.1X10"® 

-.18x10"® 

20.6x10"® 

37. lxio"® 

Modmor-I graphite/polyimide 

33. 2x10® 

228. 9 

.9- 10® 

6.2 

. 6x10® 

4. 1 

.24 

-.3X10"® 

-.54x10"® 
1 

20.5x10"® 

36. 9x10"® 


TABLE III. - THEORETICAL UNIAXIAL STRENGTHS OF VARIOUS COMPOSITES AT ROOM TEMPERATURE 
[Normalized at 0. 55 fiber volume ratio and zero voids (refs. 12 and 13). ] 


Composite system 



Uniaxial strength at room temperature 



S 

Z11T 

s mc 

S Z22T 

s 

Z22C 

S U2S 

ksi 

MN/m 2 

ksi 

MN/m 2 

ksi 

MN/m 2 

ksi 

MN/m 2 

ksi 

MN/m 2 

Boron/aluminum 

176 

1214 

185 

1276 

14.7 

101 

16.1 

111 

15.9 

110 

Boron/epoxy 

214 

1476 

188 

1296 

8. 6 

59.3 

30.1 

208 

13.6 

93.8 

S-glass/epoxy 

204 

1407 

143 

986 

4.4 

30.3 

28.3 

195 

8.3 

57.2 

Thornel-50 graphite /epoxy 

106 

731 1 

56.6 

390 

3.7 

25.5 

17.6 

121 

2.6 

17.9 

Modmor-I graphite/polyimide 

138 

j 952 

136 

938 

6. 1 

42. 1 

19.7 

136 

7.8 

53.8 




CS- 52405 


Typical fiber composite geometry (plane ABCD is the reference plane) 




Fiber volume 

ratio 0-35 



30 60 90 0 30 60 


Ply orientation angle, 0, deg 

(a) Stress in the 0 C ply. (b) Stress in the +0 ply. 



Ply orientation angle, 0, deg 


(c) Stress in the +0 ply. 

Ply residual stress for boron/aluminum composite. Composite geometry, 8[2(0), 2{±0), 2(+0), 2(0>1 ; temperature difference, 900° F (773 K) 



Ply 

orientation, 



0 30 60 90 


Ply orientation angle, 6, deg 

Figure 4. - Ply residual transverse stress for boron aluminum composites. Com- 
posite geometry, 8[0, 2(+6), 2(90), 2(^6), 0] ; fiber volume ratio, 0.55; tempera- 
ture difference, 900° F (773 K). 



Ply orientation angle, 0, deg 


Figure 5. - Outer-ply residual transverse stress for boron aluminum composites. 
Fiber volume ratio, 0.55; temperature difference, 900° F (773 K). 


I 


20 


Transverse stress, ksi 



55 ^ // — ■' " 

/ / -.75 
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0 30 60 90 

Ply orientation angle, 0, deg 

(a) Stress in the ±6 ply of composite 8[2(0), 2(±0), 2(+0), 2(0)]. (b) Composite geometry, 8[0, 2{±0), 2(90), 2(^0), 0] ; Fiber volume 

ratio, 0,55. 
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Transverse stress, ksi 
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Fiber volume 
ratio 

0. 35^ 






Ply orientation angle, 0, deg 


(a) Stress in the +0 ply of composite 8[30), 2(±0), 2(^0), 2(0)]. (b) Composite geometry 8[0, 2l±0), 2(90), 2(T0), 0] 



0 30 60 90 


Ply orientation angle, 0, deg 


(c) Outer ply; fiber volume ratio, 0.55. 

Figure 7. - Ply residual transverse stress for S -glass/epoxy composites. Temperature difference, 300° F (440 K) 
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Ply orientation angle, 0, deg 


(a) Stress in the +6 ply of composite 8 [2(01, 2(+0), 2(T0) f 2(0)]. 
10 — 


(b) Composite geometry, 8[0, 2(±0), 2(90), 2pr0), 0]; fiber volume 
ratio, 0.55. 
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Ply orientation angle, 0, deg 


co 

03 


(c) Outer ply; fiber volume ratio, 0.55. 

Figure 8. - Ply residual transverse stress for graphite Thornel -50/epoxy composites. Temperature difference, 300° F (440 K). 
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Fiber volume ratio, kf 

(a) Transverse stress. 
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Composite geometry 

-24 [4(430), 2(415), 12(0), 2(415), 41*30)] 
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(b) Shear stress. 
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(c) Longitudinal stress. 


Figure 9. - Ply residual stress for Modmor-I/polyimide composites. Temperature difference, 600° F (606 K). 
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Figure 10. - Effects of voids on ply transverse residual stress. Modmor-I/polyimide 
composites; fiber volume ratio, 0.50; temperature difference, 600° F (606 K). 



Fiber volume ratio 


Figure 11. - Effects of variable fiber volume ratio on ply residual stress. 
Modmor-I/polyimide composite; composite geometry, 24 [4 (±45), 2{±22. 5) 
12(0), 21*22.5), 4(*45)] ; void volume ratio, 0; temperature difference *600° F 
(606 K). 
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Figure 12. - Effect of having different groups of plies at different fiber 
contents on ply transverse residual stress. Modmor-I/polyimide 
composite; composite geometry, 24 [4(±45), 2(±22.5), 12(0), 2(+22.5), 
4(^45)]; zero void content; temperature difference, 600° F (606 K). 
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Matrix thermal coefficient of expansion, in./ in. /°F 
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Matrix thermal coefficient of expansion, cm/cm/K 

Figure 13. - Effects of matrix thermal coefficient of expansion on ply transverse 
residual stress. Modmor-I/polyimide composites; zero void content; fiber 
volume ratio, 0.50; temperature difference, 600° F (606 K). 
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Figure 14. - Effect of matrix modulus on ply transverse residual stress. Modmor-I/ 
polyimide composites; zero void content; fiber volume ratio, 0.50; temperature 
difference, 600° F (606 K). 
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Figure 15. * Adjacent ply relative rotation dueto residual stress for 
graphite Modmor-I/polyimide composites. Temperature difference, 
600° F (606 K). 
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Composite geometry 

24[6(±45), 12(0), 6(^45)] 
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Figure 16. - Transitional plies and decreasing number of off-axis 
plies effect on ±45° ply residual transverse stress. Modmor-I/ 
polyimide composite; temperature difference, 600° F (606 K). 
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